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ABSTRACT 

The ability of the neutron star surface to supply all or only part of the charges 
filling the pulsar magnetosphere is crucial for the physics prevailing within it, with 
direct consequences for the possible formation of pair creation regions. We evaluate the 
conditions for e~ emission from pulsar surfaces for a simple Goldreich- Julian geometry 
taking both thermal and field emission processes into account. Using recently published 
estimates for the equation of state at the neutron star's surface, we show, that for a 
large range of T sur f , B and P, the liberated charges will fully screen the accelerating 
B-parallel electric field En. For surface temperatures T sur f < 2 • 10 5 K a balance between 
field emission of electrons and shielding of the field will occur. Even in the overidealised 
case of T sur f = one can expect a prodigious supply of electrons which will weaken 
the accelerating Em. We calculated the motion of electrons along selected polar field 
lines numerically for the low temperature, field emission scenario yielding their Lorentz 
factors as well as the produced radiation densities. Inverse Compton and later curvature 
losses are seen to balance the acceleration by the residual electric fields. We found that 
the conditions for magnetic pair production are not met anywhere along the field lines 
up to a height of 1500 pulsar radii. We did not a priori assume an "inner gap", and 
our calculations did not indicate the formation of one under realistic physical conditions 
without the introduction of further assumptions. 
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1. Introduction 

A proper assessment of the conditions very close 
to the surface of a neutron star is a vital ingredient 
for any model that attempts to explain the enigmatic 
origin of the observed electromagnetic radiation from 
pulsars. 

Numerous studies of the polar cap acceleration 
have been made within a particular framework: A 
steady acceleration of electrons within a space charge 
limited particle flow originating at the surface of the 
neutron star will produce particles emitting magnet- 
ically convertible curvature 7— rays (e.g. Arons & 
Scharlemann 1979). These high energy photons are 
converted into e + e~-pairs with high enough efficiency 
(Daughcrty & Harding 1982) so that in the end, they 
and their daughter products will exceed the initial 
(Goldreich- Julian density) by a factor of 10 4 (e.g. 
Zhang & Harding 2000). An important a priori as- 
sumption is inherent in all these calculations: Above 
a critical height /i gap a region of high conductivity 
will quench any electric field Ey parallel to the guid- 
ing B-ficld lines, so that the full potential across the 
polar cap will be active between the outer surface of 
the neutron star and the inner surface of the high 
conductivity region. 

Our approach differs in that respect from other 
workers, we do not make that initial assumption, but 
try to calculate the amount of charges produced by 
the surface, how they affect and are affected by Ey, 
their energy losses and the possibility of pair produc- 
tion. Should one find suitable conditions for such a 
pair production region to exist, the methods quoted 
above will be the prime choice for its description. 

The conclusion that electron-positron pairs with 
high densities exist close to the neutron star has re- 
cently been questioned by several authors. Kunzl 
et al. (1998) considered this widespread assumption 
of dense pair plasmas for the problem of radio wave 
propagation within the neutron star magnetosphere. 
On hand of the most recent emission height estimates 
(50 - 100 i?Ns) for the radio emission (Kijak & Gil 
(1997); Kramer et al. (1997, 1998)) they found that 
the existence of low frequency radio emission for ex- 
ample from the Crab pulsar is incompatible with the 
scenario of extended dense pair production. These 
findings were fully confirmed by the theoretical in- 
vestigation of Melrose & Gedalin (1999), who treated 
the plasma emission mechanisms in detail and came 
to the conclusion that only low energy particles with 



Goldreich- Julian density could exist in the radio emis- 
sion regions in order to alow for escaping radio waves 
(see also Melrose (2000) and references therein). Fur- 
thermore recent simulations of the pair cascade by 
Arendt & Eilek (2000) for strong and weak magnetic 
fields shed some light on the difficulties for an effi- 
cient production of secondary pairs by high energetic 
primary particles via 7-ray curvature photons. Al- 
though their findings need some more investigation 
they could show that the field strength at the sur- 
face has the most important effect on the cascade. 
For strong fields B ~ 10 9 T they find almost no sec- 
ondary photons and the cascade terminates in one 
generation. The weak fields (B ~ 10 7 T) had either 
no cascades or quite weak ones. Only for moderate 
fields (B <~ 10 8 T) they observed some secondary par- 
ticle generations. Arendt & Eilek concluded that it 
is not yet clear whether pulsars can really produce a 
large number of pairs at the surface gaps. 

This tempted us to question the canonical a priori 
assumption of the existence of an inner pair creating 
gap and investigate its possibility given a reasonable 
estimate of the surface characteristics and the fields 
involved. 

In the literature there are three fundamentally dif- 
ferent ideas for pair production in the inner magne- 
tosphere: A "starved" magnetosphere where the neu- 
tron star surface provides only an insufficient supply 
of charged particles (e.g. Ruderman and Sutherland 
(1975)), a magnetosphere where the Goldreich- Julian 
density is reached at the surface but due to curvature 
effects of the field lines deviations occur close to the 
surface (e.g. Arons & Scharlemann (1979), Barnard 
& Arons (1982)) and Shibata's model which predicts 
electric fields caused by overdense regions, that means 
in the opposite sense than in the Arons model (Shi- 
bata, 1997). 

The latter two ideas suffer from serious weaknesses 
since it turned out that Arons type models do not 
work for slow pulsars (Beskin (2000) gives an upper 
limit of 0.1 to 0.3 s for the period). The Shibata model 
in turn neglects the back reaction of an overdense 
region which reduces the particle flow from inside. 
Therefore we will restrict our discussion to models 
of Rudcrman-Sutherland type and fix the parameter 
range in which this kind of models is reasonable. 

For the scope of this article we consider a strong 
magnetic dipole anchored in a heavy (1.4M sun ) sphere 
with a radius r = 10 km that rotates twice per second 
(Q = 2ir ■ 2s~ 1 ). Let the B-field have a strength of 
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I .Bo I = 10 8 T at the surface and for simplicity's sake 
let the magnetic and the rotational axes be aligned^. 



B(r,0) = B, 



(cos 9, - sin 0, 0) (1) 



Such objects have been discussed as idealised cases of 
radio pulsars by Goldreich & Julian (1969) and many 
authors since then. Of course one has to realise that 
the above mentioned idealisations limit the applica- 
bility of such simple models for the interpretation of 
observational facts from the growing number of pul- 
sars that is being discovered. We will however try to 
extend the scope of these models by constraining the 
number of charged particles available and analysing 
the fate of individual charged particles (limiting our- 
selves to electrons and positrons) in the polar regions 
of such heavy rotating dipoles under the constraints 
given by some of the more recent observations. Al- 
though a pulsar magnetosphere is characterized by 
extreme physical conditions in terms of strong elec- 
tric and magnetic fields the fundamental processes 
that can supply charges are reasonably well known 
from terrestial experiments and technical applications 
and can be extrapolated to the surface conditions of 
a neutron star. Such an extrapolation provides con- 
straints on the initial charge density which is vital 
for any model of cascading particle creation and the 
radiation processes. 

2. Accelerating fields and charge densities 

Muslimov & Tsygan (1990, 1992) showed how the 
electric and magnetic fields of a rotating neutron star 
are distorted by relativistic effects. A comprehen- 
sive treatment of the rotating dipole in the frame of 
general relativity was given by Muslimov & Harding 
(1997). For our purposes of demonstrating the ma- 
jor physical properties near the polar cap it suffices 
however to stick to a classical description (including 
special relativity of course). 

Our rotating magnetic dipole acts as an unipolar 
generator and induces an electric field as the conse- 
quence of its surface charge distribution 
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(3cos 2 0- 1, 2 cos sin 0,0). (2) 



1 We use spherical coordinates (r,8,<j>) centered in the neutron 
star with 9 being the directional angle referred to the pole which 
coincides with the rotation axis. 



The component parallel to the magnetic field lines 
E ll = W is given by 
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and reaches 2.5 • 10 13 V/m at the surface near the 
magnetic poles (Goldreich & Julian 1969; Holloway 
1973). In vacuum all appropriately charged particles 
will be quickly accelerated to very high energies, but 
because of the strong magnetic field they are con- 
strained to move along the magnetic field lines given 
by the parametric equation 



r(0, O ) = r NS 



sin 2 
sin 2 0q 



(4) 



Here 0o is meant to be the colateral angle of the field 
line at the surface of the sphere. 

Goldreich & Julian (1969) already calculated the 
electron density which shields the electric field and 
yields the only stable static solution for a pulsar mag- 
netosphere (excluding centrifugal and gravitational 
forces) . 
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with e being the electron charge. 

In our example this amounts to tt-gjo = 2eo ^ B ° = 
1.4 • 10 17 m~ 3 at the poles on the surface. Deutsch 
(1955) was the first who showed that the high induced 
fields make it unlikely that stellar rotators with strong 
magnetic fields will rotate in a (charge-free) vacuum. 
Several well studied processes enable i.e. electrons to 
escape from a metallic surface. If the neutron star 
surface can thus liberate sufficient charges to keep the 
Goldreich- Julian current flowing, the emission will be 
space-charge limited like an electronic valve (diode). 
Electrons will be accelerated awayQ from the polar 
regions if Ci ■ B > 0F|. 



2 The observed high brightness temperatures of the pulsar radio 
emission suggest the operation of beaming mechanisms for the 
enhancement of the radio intensity. We'll therefore focus our 
attention on outflowing particles where intensity enhancement 
by forward beaming would be natural. 

3 In some models in the literature we need anti-alignment of 
Q and B. In that case the Goldreich-Julian charge density 
above the polar cap is positive. We do not discuss this case 
here but since the binding energies of ions are assumed to be 
comparable to those of electrons the qualitative argumentation 
remains unchanged. 



3 



Usov & Melrose (1995) have reminded us of the 
fact that the relevant fundamental surface emission 
processes are field emission and thermal emission 
(perhaps enhanced by the Schottky effect (Schottky 
1914)). Thermal surface emission of electrons was 
discovered by T.A. Edison (1884) and has been ex- 
tensively studied since the beginning of our century 
(Dushman 1930). Field emission from metal surfaces 
is also known for a long time (Fowler & Northeim 
1928) and used in many technical applications. The 
effectiveness of both emission processes is controlled 
by £\y, the work function or binding energy, the en- 
ergy required to liberate an electron from the crys- 
tal lattice at the surface boundary. Values for £\y 
and the Fermi-energy Ep for the pulsar surface can 
only be extrapolated from terrestial values as the high 
densities of pp c ~ 10 6 kg m~ 3 and magnetic fields 
B ?s 10 8 T cannot be reached in laboratory experi- 
ments. We will therefore first try to establish a rea- 
sonable range of values for E^ in order to investigate 
the effectiveness of surface emission. 

2.1. Constraints for the range of the work 
function 

Using a Thomas-Fermi-Dirac-Weizsacker approxi- 
mation Abrahams & Shapiro (1992) calculated the 
cold (T = 0) surface density pp c to be about 2.9 • 
10 6 kgm~ 3 . That result enables us to give an im- 
proved estimate of the Fermi energy at the surface: 



Ef(pfc) = 



2-tt 4 H 4 c 2 /ppc-(28-2) 



e 2 B 2 m e 



56 • m c 



(6) 



which amounts to Ep = 4.17 • 10 2 eV (as usual Ti = 
/i/2"7r). In our subsequent calculations we take that as 
the most likely value of Ep . In the case of an unknown 
work function -Ew, one usually approximates £\y = 
Ep. So far calculations of the actual electron work 
function have not been made and for want of a better 
alternative we will use the Fermi energy as the best 
available measure of the work function. Abrahams & 
Shapiro (1992) give the ionisation energy for iron as 
120 eV and we will consider the range from ionisation 
energy to cyclotron energy (0.12 - 12 keV) to be also 
a reasonable constraint for £\y- 

2.2. Field (cold cathode) emission 

By calculating the transmission coefficient of elec- 
tron wave functions through a potential wall 



U(x) = e- 1 E w - E • x for x > and U(x) = for 
x < (E being the electric field) Fowler & Nord- 
heim (1928) provided a description of electron cur- 
rents from metal surfaces in good agreement with the 
experimental results and useful in many technical ap- 
plications since then. Using their expression of the 
cold cathode current density and assuming a relativis- 
tic flow of electrons we find for the particle density 



"field (Eg) = 



/Ew . ip2 
E F ^0 



2-kHc (E w + Ep) 



In Fig. 1 we plot the ratio of the field emission 
charge density to the Goldreich-Julian density for a 
standard pulsar. The authors like to emphasize the 
point that this density is independent of the some- 
what uncertain surface temperature and depends only 
on Eq. 
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Fig. 1. — Cold cathode emission n/riQj versus work 
function E-w for a standard pulsar (see text). 

Field emission can supply the Goldreich-Julian 
density for a cold neutron star surface for any Ew < 
2keV! Fig. 2 however shows, that for the extreme 
case of a 100 times weaker field (i.e. a slowly rotat- 
ing weak pulsar) the maximum binding energy drops 
down as far as 50 eV. 

Fig. 3 shows how for an assumed work function of 
i.e. 500 eV, the cathode efficiency increases dramati- 
cally with the electric field. We may therefore expect 
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Fig. 3. — The cold cathode efficiency n/ncj at E^i 
417eV as a function of the surface field E . 



Fig. 2. — The critical work function for the Goldreich- 
Julian density versus the surface field E . 

field emission to play a significant role in the supply 
of charges for the average pulsar magnctosphere. 

2.3. The thermal electron emission 

In contrast to the quantum mechanical effect of 
field emission, the thermal surface emission of elec- 
trons from a conductor is a purely classical effect: For 
a given particle distribution function, be it Maxwellian 
or Fermi-Dirac, there will be a fraction of particles 
with a kinetic energy that exceeds that of the po- 
tential barrier characterised by the binding energy. 
That number of electrons liberated from a hot metal- 
lic surface is given by Dushman's equation (Dush- 
man 1930). In the presence of strong external electric 
fields, the number of electrons is increased over that 
given by Dushman's equation because of the Schottky 
effect (Schottky 1914) which lowers the effective po- 
tential barrier for the particles. The charge density is 
thus described by the combined Dushman-Schottky- 
Equation 

n BS (T,E ) = ^-T 2 -c - + -V 4 ^ (8) 
e • c 

(k being the Boltzmann constant) with A a = 47r -"^ e - fc2 
= 120.2 Acm -2 which depends only on £\y In Fig. 
4 we show the contributions to the current density 



according from field emission and from thermal emis- 
sion with and without the Schottky effect for a range 
of temperatures likely for a standard pulsar with an 
assumed E^j = 417 eV. 

As has been pointed out before, for E n = 2.5 • 10 13 V/m 
field emission alone can supply all that is required, 
but for temperatures above 10 5 K, thermal emission, 
which is only weakly dependent on the effective field 
becomes even more efficient. From 2 • 10 5 K onwards, 
thermal electrons, emitted from the surface even with- 
out a large electric field can supply the GJ-density. 
Fig. 5 brings everything into perspective by outlining 
the extremes of the parameter space. 

We show the effective ratio of total charge emission 
for work functions of Ew = 2065eV and £\v = 417 eV 
as well as for a field of y^. Only for extremely high 
E w > 1200 eV or very low fields E < 10 12 V/m and 
low surface temperatures T < 2 • 10 5 K the electron 
supply from the pulsar surface will be insufficient to 
fill the magnctosphere with charges to shield the in- 
duced electric field. It is however doubtful that ob- 
jects with such low temperatures can be found as in- 
ternal stresses caused by the rotational braking of 
the neutron star set a lower limit of a few 10 5 K 
for the surface temperature (Tsuruta 1998 and refer- 
ences therein). But normally the potential supply of 
charges exceeds the requirements of the G-J-density 
by a wide margin. Thus any accelerating electric field 
at the surface of the neutron star will be quenched 
within an extended corotating charged lower region of 
the magnetosphere. Therefore we do not expect any 
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Fig. 4. — The surface efficiency n/riQj at £\y = 
417cV for a standard pulsar as a function of the sur- 
face temperature. Thin solid line: thermal emission, 
dots: thermal emission with Schottky effect, Thick 
solid: cold cathode emission. 



significant acceleration in the inner magnctosphere 
which, of course, implies the absence of highly rel- 
ativistic particles and no pair production at all under 
these circumstances. 

3. Neutron stars dominated by field emission 

We have shown in the last section, that high bind- 
ing energies or low surface temperatures are the only 
conditions compatible with significant surface electric 
field strengths that can accelerate particles. In these 
cases, an equilibrium field of E cq = E • (1 — will 
remain when n e electrons flow off the surface. Both 
the electric field Eh and the shielding density are pro- 
portional to the product B • Q. We introduce a scale 
factor £ = rF jp c " n to account for the possible empiri- 
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cal variation of B • Q and by numerically solving the 
balance equation 
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Fig. 5. — The surface efficiency n/riQj for different 
Eyj and as a function of the surface temperature, long 
dash: total efficiency for Ew = 417eV, solid: total ef- 
ficiency for Ew = 417eV and O.lEo, short dash: total 
efficiency for Ew = 2065eV, dots: total efficiency for 
Ew — 417eV and a weak field pulsar with 0.01E and 
0.01n GJ . 

for 0.02 < £ < 5 we obtain the equilibrium densi- 
ties for any given temperature T and work function 
Ew covering a wide range of neutron stars with dif- 
ferent B and f2. Fig. 6 shows the variation of the 
residual electric field £E(1 — -jr^) as a function of 
£ for E w = 417cV and E w = 1200eV. A tempera- 
ture of 1.85 • 10 5 K, just below the onset of thermal 
emission for Ew — 417eV, was chosen for the calcu- 
lation. The residual field grows with £ until it satu- 
rates at 2.1 • 10 12 V/m for £ = 1. A higher Ew causes 
earlier saturation at a higher level. Between 300 eV 
and 2000 eV the residual field saturation level E\i m 
depends nearly linearly on Ew and can be approxi- 
mated as E Um w 1.3 ■ 10 10 • {E w - 250eV)V/m. For 
our standard Ew, the density n deviates by only 10% 
from the Goldrcich- Julian value, but higher Ew show 
lower initial densities with the expected increase in 
charge density with increasing electric fields. 

3.1. The equilibrium energy of electrons in 
the magnetosphere 

Any residual field will quickly accelerate the elec- 
trons to relativistic energies, balanced by the braking 
due to various energy loss processes. Because of the 
high magnetic field strength, the motion of charges 
will be along the field lines. With the help of the 
parametrisation (4) we can describe a particle tra- 
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6. — The residual electric surface field 
) at T = 1.85 • 10 5 K as a function of £ for 



jectory using only its co-lateral angle at the pulsar 
surface (starting angle 6 S ) and the co-lateral angle 9 
on the trajectory. At each point of such a trajectory, 
the energy gain of a charge by the parallel component 
of the residual electric field 
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m e c 



(1 



)E||(Mo) (10) 



will be balanced by a number of energy loss mecha- 
nisms. Curvature radiation is an obvious loss mecha- 
nism for relativistic charged particles forced to follow 
a curved trajectory. The energy loss of a relativistic 
particle due to curvature radiation is given by (e.g. 
Zheleznyakov 1996, p.231) 



dj 2 4 e 2 1 

dt ' CR ~ 3 7 4ire Q m e cR c (9,9 Q ) 2 



(11) 



Here we use the curvature radius of a dipolar field line 
in its parameterised form (Lesch et al. 1998) 



R c (9, 9q 



r NS sin0 (l + 3cos 2 fl) 
K 1 + cos 2 9 



3/2 



sin tf 



(12) 



The curvature radii in the polar regions of a neutron 
star with a dipolar magnetic field are of the order of 



10 6 m. To balance the huge acceleration of an un- 

1.5 • lO 1 ^- 1 by 
l 



shielded field at the surface ^ 

dt 



8.5 • 10~ 19 7 4 s" 



curvature radiation with |cr 
one needs Lorentz factors of 7 w 3.6 • 10 8 . Particles 
with smaller Lorentz factors experience only insignifi- 
cant losses due to the 7 -4 dependence of the curvature 
losses. 

Inverse Compton scattering of relativistic particles 
on the thermal photons emitted by the hot neutron 
star surface (Sturncr 1995) limits the electrical accel- 
eration of particles near the pulsar surface (Supper & 
Triimper, 2000). If one also includes the interaction 
of the charged particles with their own production of 
scattered photons one would expect an even more ef- 
ficient braking to take place. Due to the curvature 
of their trajectories and we expect a fairly strong in- 
teraction of the particles with their own previously 
upscattered radiation within a few additional meters 
of travel, so that the effective energy density of the ra- 
diation field will increase quickly in the particle frame. 

For an assumed complete dissipation of the avail- 
able energy of the residual electric field near the sur- 
face we can estimate a effective temperature of 



T = 



= 7.3-10 7 K(£ )5af 

™GJ 



(13) 



as a natural limit for the produced radiation. Here we 
use the blackbody constant a = 7.566- 10 _16 Jm~ 3 K~ 4 

and x — — ^— as the normalised distance from the sur- 

. rNS 

face. It is however unrealistic, to expect such efficient 
conversion of the available energy into radiation. Fol- 
lowing Sturner (1995) we expect dissipation scales for 
the dominant processes of the order of a few m in the 
vicinity of the neutron star surface. Again a simple 
estimate from a balance of the energy gain due to the 
electrical field on a dissipation scale A will give us an 
improved temperature limit 



eE|,(l )n\ 



4-10 6 Ka;- 1 (1 )n\ 



(14) 



which turns out to be compatible with the observed 
x-ray temperatures of pulsars (Tsuruta 1998). 

A more reliable description of the conditions near 
the pulsar surface can be attained by computing the 



7 



Lorcntz factor of electrons with 
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(15) 



together with the development of the energy density 
of the radiation field 



along a given field line parameterised by the colat- 
eral angles 9 and 9 from the surface out to 10 3 rNs- 
The energy losses by the major processes outlined 
by Sturner (1995) as there are: curvature radiation 
(^g Icr ), the Klein-Nishina correction as well as mag- 
netic, resonant and non-resonant inverse Compton 
scattering, were included in the term d7 ^° ff " > linvCompt • 
The self-enhanced scattering from the radiation pro- 
duced in lower levels of the magnetosphcre was ap- 
proximated by a local calculation of the effective tem- 
perature of the radiation bath T e g = (^f) 1 ^- We did 
not calculate the exact spectrum of the radiation bath 
for reasons of efficiency. But we expect our crude ap- 
proximation to be valid due to the curved trajectories 
of the electrons and the straight, tangential propaga- 
tion of the emitted photons. After a short propaga- 
tion, the scattering will not be back on, but at an 
slowly increasing angle. We repeated the calculations 
also without this enhancement effect, only to find, 
that the initial energy plateau did not extend quite 
as high up into the magnctosphere and the maximum 
energy achievable by the particles was a little higher. 
Because of the widely varying scales of the problem, 
the integration was performed for twelve segments of 
the trajectory using an Bulirsch-Stoer algorithm with 
decreasing resolution from one segment to the next 
(Press et al. 1992). The final values of each segment 
served as initial conditions for the subsequent ones. 
The results of each segment were carefully checked 
for their resolution independence, too coarse resolu- 
tions showing up as cither oscillatory behaviour or 
non-convergence. Because of the enormous accelera- 
tion of a particle starting with 7 = 1 at the surface, 
the first (surface-) segment had to be as short as 10 
cm subdivided into 200 steps. Most interesting was 
the limiting case of T = 2- 10 5 K and n = 0.9n GJ . Fig. 
7 provides a plot of the achieved Lorentz factors and 
effective temperatures as a function of distance from 
the surface. After travelling about 5 cm, the acceler- 
ation was balanced by non-resonant inverse Compton 
losses at 7 — 2.2 -10 5 . After that initial peaking of the 



Lorentz factors, the particle energies decreased slowly 
down to 7 = 1.1 • 10 5 at h = 20m acccompanied by a 
slow rise of the effective temperature. That decrease 
is a consequence of self-enhancement and does not 
happen without it. Over the next 630m the effective 
temperature rises quickly to 10 6 K and the Lorentz- 
factor reaches 1.3 • 10 6 . Further, slower progress leads 
to a complete balance and saturation at h = 1.3 - 10 5 m 
with T = 2.9 • 10 6 K and 7 = 1.1 • 10 8 . At 5.5 • 10 4 m 
curvature radiation takes over as the dominant en- 
ergy loss mechanism. The Lorentz factors now de- 
crease slowly down to 7 = 3.9 • 10 7 at a distance of 
1.5 • 10 7 m. 
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Fig. 7. — The Lorentz factors (thin solid line: with 
self-enhanced inverse Compton scattering, thick solid: 
without enhanced scattering) and effective tempera- 
tures (dashed line) for enhanced scattering as a func- 
tion of distance from the surface for T = 2 • 10 6 K and 
n = 0.9 • nGj ■ 

A number of different calculations with varying ini- 
tial parameters n < 0.1 • uqj and T were made, but 
showed no qualitative difference in the behaviour of 
their solutions. The case described can hence be taken 
as a reasonable representative of a " starved" system, 
having a low surface temperature and producing a 
strong source of X- and 7- rays with T e g w 3 • 10 6 K. 

3.2. Pair production 

Our calculations differ from those of Sturner (1995) 
in as far as we did not a priori assume the existence of 
an inner polar gap. It is often unquestionly assumed 
that after about 100m a particle creation avalanche 
creates a dense n 3> uqj charge neutral plasma that 
shorts the available electric field. For an estimate of 
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the location of such an avalanche region we calculated 
the magnetic pair creation probability for curvature 
photons of a characteristic energy 



E F . 



3chj 3 



(17) 



and Compton boosted thermal radiation having an 
energy of 

Bcomp = l 2 kT c s (18) 

The hard photons were assumed to be emitted tan- 
gentially and travelling outward from the field lines 
over a distance A, their angle with the magnetic field 
lines 0(A, 0, 0q) and the corresponding component of 
B perpendicular to the path of the photon 

B± (A, 0, 6 ) = B ■ x~ 3 sin 6(A, 0, O ) (19) 

were explicitely calculated at each step along the tra- 
jectory of the accelerated charge. We used the local 
strength parameter 



x( 7 ,A,Mo) = 



2 Bcrit 



(20) 



and Erber's (1966) Besselfunction approximation^] to 
calculate the local attenuation factor 



a 7 B 



a 0.16 ( 2 



(21) 



with a = the electron Compton wavelength 

A e = -2- and the critical magnetic field £> cr it = 

2^! = 4.414 . iq9t. K 1/3 stands for the modified 
Besselfunction of the third kind and order 1/3. Mul- 
tiplying the attenuation factor a 7 s with the local 
steplength of the integration provides the local pair 
creation probability p 7 b(A, #o) which in turn mul- 
tiplied with the number of emitted photons, ie. 

2e 2 



9cheoR c 



■7 = 2.1 • 10" 8 7m" 



(22) 



for curvature radiation, gives the local number of 
pairs created per accelerated chargefl. We tested the 



4 Although that approximation involves more effort for its eval- 
uation as compared to the commonly used exponential form, it 
accurate over a greater range of \ and, having no discontinu- 
ities, more useful for numerical work. 

5 Even under favourable conditions there will be only a small 
number of of initial pairs produced. This is a consequence of 
the fact that p 7 s < 1 and that, at least for curvature radiation, 
the number of productive photons per unit length is very small. 
Because of the lower secondary and tertiary Lorentz factors, the 
conditions for efficient pair production in the form of daughter 
products tend to be even worse. 



algorithm for the physically unrealistic case of miss- 
ing inverse Compton braking and did indeed find a 
pair production zone at a height of 200 m above the 
pulsar surface. 

However, of all our simulations that included non- 
resonant inverse Compton scattering, none showed 
PjB > anywhere near the pulsar surface or within 
the simulated range of distances, neither for curvature 
photons nor for maximally boosted inverse Compton 
photons. This is due to the fact that the calculated 
Lorentz factors do not reach their highest values near 
the pulsar surface but further out, where the magnetic 
field has already fallen drastically due to its a; -3 de- 
pendence. 

4. Discussion and conclusions 

We have seen that the surface temperature and the 
work function play a crucial role in determining the 
initial energy and density of the lower pulsar magne- 
tosphere. Two scenarios are possible: Copious ther- 
mal emission of electrons leading to a "saturated" 
charged magnetosphere with n = uqj and E|| = 0. 
In that case we expect the inevitable potential drop 
and energy dissipation further away from the pulsar 
surface. Alternatively a "starved" magnetosphere, 
where n < riQj is determined by the residual equi- 
librium field (1 — )En. Here our simulations differ 

V "GJ ' , U 

from those made previously by others as we include 
the shielding of most of the E\\ by the field emission 
produced particles together with the inverse Comp- 
ton loss processes. By taking both physical effects 
into account, we find that conditions for pair creation 
cannot be met near the surface of ordinary rotating 
neutron stars that form most of the observed pulsars. 

In the first case, we do not expect any Lorentz 
factors 7 > 1. We have, however only described 
an unperturbed equilibrium atmosphere. The exten- 
sively studied fact that pulsar radio emission varies 
on timescales down to microseconds (Lange et al. 
(1998)) can be seen as an indication that such an 
atmosphere has inherent instabilities. Theoretical ex- 
amples of instabilities in saturated magnetospheres 
were first given by Mestel et. al. (1985). Because of 
the strength of the compensated electric field, even 
smallest perturbations in the surface densities will be 
accompanied by (mildly) relativistic particles. Kunzl 
et. al. (1999) have described a scenario in which such 
perturbations lead to the radio emission via a drift 
induced free electron maser process. We expect that 
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most neutron stars to fall into that category, having 
surface temperatures higher that 2 • 10 5 K combined 
with an E\y ~ 400cV. If, as proposed by Pavlov 
& Zavlin (1998) neutron stars are surrounded by a 
gaseous atmosphere, we may expect easy liberation 
of charges in the form of electrons and ions, and a 
saturated magnetosphere is then unavoidable. 

In the second category we would find cold neutron 
stars with a relativistic X-ray emitting lower magne- 
tosphere, having effective temperatures of a few 10 6 K. 
The achieved Lorentz factors are of the order of sev- 
eral 10 5 in the lower magnetosphere reaching a peak 
of 7 = 10 8 at about 10r N s and decrease monotoni- 
cally to a few 10 7 at 1500r NS . A significant flow of X- 
and 7— rays can be expected in these circumstances, 
but the authors like to emphasise that radio emission 
models generally require much lower Lorentz factors 
(w 10 — 1000). The discussed loss processes are evi- 
dently incapable of braking the charges hard enough 
in order to achieve this, making difficult to explain 
any significant radio emission in these circumstances. 

But the Lorentz factors in the region of strong mag- 
netic fields are nevertheless too low to expect signif- 
icant pair production and the subsequent formation 
of an inner gap in these cases. 

But pulsars are radio and x-ray sources at the same 
time, this fact can only be reconciled with the ar- 
guments above by fully accepting the observational 
evidence that the magnetosphere is far from a homo- 
geneous or steady state. A global investigation of in- 
stabilities in the current system of neutron star, mag- 
netosphere and its outer boundary ought to be un- 
dertaken. Ideally, it would show that spatial or tem- 
poral variations in the magnetosphere would create 
the observed x-rays and the radio emission in different 
parts of the perturbation. Our calculations are meant 
to cover the possible extremes of the surface emis- 
sion and the resultant charge densities and Lorentz- 
factors. Our terrestial experience show us, that sur- 
face electron emission is affected by random variations 
of emission efficiency (formation of emitting islands) 
that are mainly responsible for the low-frequency part 
of a diode noise spectrum. It is not unthinkable that 
the thermal emission efficiency fluctuates over the po- 
lar cap. In that case we can expect local regions of 
the lower magnetosphere also to fluctuate individually 
between the two extreme cases described. 

A few self critical remarks ought to be made: This 
study presents a drastically simplified scenario. We 
have omitted the effects of the pulsar's inclination 



of the magnetic axis and possible surface inhomo- 
geneities. Hence a greater spatial variation between 
the two extremes of starved and saturated charge lib- 
eration would be expected for real pulsars. We omit- 
ted a treatment that includes the effects of general 
relativity as i.e. was done by Muslimov and Harding 
(1997). If these effects are taken into account, fields 
and consequently Goldreich- Julian densities will be 
higher. 

But due to the steep (exponential) characteristic 
of the surface emission processes, one expects a shift 
of the critical temperature by perhaps 10 5 K. We also 
neglected resonances of the pair production process 
described by Daughcrty & Harding, (1996) which lead 
to an enhancement of the pair production mechanism. 
The Lorentz factors compatible with inverse Comp- 
ton braking however are by a factor of 10 3 lower 
than those constrained only by curvature radiation 
when pair production would be possible near the polar 
cap, and complete description of the problem would 
of course include both general relativistic effects and 
these resonances. A more detailed investigation of the 
physics of the upper surface layers would provide us 
with a more precise determination of Ew without the 
approximation of Ew = Ep. Furthermore, the sim- 
ple approximation for T e g does not incorporate the 
detailed spectral characteristics of the produced radi- 
ation, but the produced hard radiation will interact 
strongly with the accelerated electrons, we may even 
have underestimated the scattering losses. 

We are lead to conclude, that the fact that we ob- 
serve pulsars with X-ray temperatures of several 10 6 K 
excludes the formation of an inner gap, and the fact 
that radio emission is also observed from such sources 
provides another argument against the high n 3> wqj 
densities that are commonly expected from an inner 
avalanche region (Kunzl et al. 1998). We are also led 
to conclude, that observed pulsars are dominated by 
free charge emission with n = jt-qj, thus it might be 
advantageous to study the possible mechanisms for 
radio- and X-ray emission in that context. 
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